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Abstract 

Thermodynamic Bethe ansatz equations for XXZ model at |A| > 1 is simplified to an integral 
equation which has one unknown function. This equation is analytically continued to |A| < 1. 

1. Introduction 

Thermodynamic Bethe ansatz equations for exactly solvable one-dimensional systems have 
many unknown functions. HJ About the XXZ model at |A| > 1, Gaudin- Takahashi equation 
contains infinite unknown functions. p], |3|] Here I can simplify this set of equations to an integral 
equation which contains only one unknown function. For Hamiltonian 

N 1 

H(J, A, h) = -Jj2 S?S? +1 + SfSf +1 + A(5f Sf +1 --)-2hY: Sf, h > 0, (1) 
i=i 4 i=i 

thermodynamic Bethe ansatz equation at temperature T is 

, , . 27rJsinh0 

In r/i (a;) = — s(x) + s * ln(l + r] 2 (x)), 

T<p 

ln?7 i (x) = s*ln(l + ^-_i(x))(l + ry i+1 (x)), j = 2,3,..., 

lim ^ = (2) 

Here we put 

A = cosh0, Q = it /(f), s(x) = - £ sech( — 2u ^ \ s* f(x) = f s(x - y)f(y)dy. (3) 

4 n=-oo * J ~Q 

The free energy per site is 

(p J-Q J-Q cosh (p — COS(0Xj 



On the contrary new equation is 

Ji / 0/ r i r 27rJsinh0 
= 2cosh(-) + J^-^cot-[ar-j/-2ijexp[ — ai(y + «)j 



0r n .i r 27rJsinh0 \ 1 dy ._. 

+ cot - [x - y + 2%\ exp — ai(y - z ) ^^tt - ' 5 
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and free energy is given by 

/ = -Tln«(0). (6) 

Contour C is an arbitrary closed loop counterclockwize around 0. 2nQ,n ^ and ±2z + 2nQ 
should be outside of this loop. This loop should not contain zeros of u(y). It is expected that 
u(y) has not zero in region \$sy\ < 1. This equation can be calculated even at imaginary <p, 
namely, |A| < 1 case. This equation converges numerically at least T/J > 0.07. The results 
coincide with those of older equations. Then this equation unifies Gaudin-Takahashi equation 
for A > 1 and Takahashi- Suzuki equation for |A| < 
2. Derivation 
We should note that if 

g(x) = s*h(x), (7) 

it stands 

g(x + i) + g(x — i) — h(x). (8) 

The Fourier transform of (0) is 

g{uj) = h(u), IjJ = ^-n. 



Then we have 



(e» + e-»)~g{u) = ~h{u), 



and obtain (H). 

Then set of equations (fj) is rewritten as 

/ -\ / -x r27rJsinh0 ^ r . /., , ^^ 

r]i{x + i)r)x(x -%) = exp — ) i i - 2nQ) (1 + 772(2;)), 

1 T(j) V J 

7]j(x + i)r}j(x - i) = (1 + 77j_i(x))(l + rj j+1 (x)), j = 2, 3, 

lim ^ = 0) 

At J = this set of equation becomes a difference equation and we have analytical solution, 

_ / sinh(j + l)h/7\ 2 

We can expand perturbatinally as power series of J/T. Very surprisingly i]j(x) has singularity 
only at x = ±ji, ±(j + 2)i + 2nQ. So we assume that f]j{x) is univalent on the complex plane 
of x and 1 + i]j(x) is factorized as follows: 

1 + rjjix) = Aj(x - ji)A~(x + j^Bjix - (j + 2)i)B~{x + (j + 2)i). (10) 

Here functions Aj(x), Bj(x) are periodic with periodicity 2Q and have singularity only at 
x = 2nQ and 

A~(x) = Aj(x), B~{x) = Bj(x). 



One should note that Aj(x) is not an analytic function of x but Aj(x) is analytic. Then second 
equation of @ becomes 

rjj(x + i)rjj(x — i) = 



Aj-i(x - (j - l)i)A,-_i(x + (j - - (j + l^B^ix + (j + l)i) 



xA j+1 (x - (j + l)i)A i+ i(a; + (j + l)i)B j+1 (x - (j + 3)i)B j+1 (x + (j + 3)i). (11) 

If we assume that rjj (x) is factorized as 

r]j(x) = X(x - ji)Y(x + ji)Z(x - (j + 2)i)W(x + (j + 2)i), 

rjj(x + i)r)j(x — i) is 

X(x - (j - l)i)Y(x + (j - l)i)X(x - (j + l)i)Y(x + (j + l)i) 
Z(x - (j + l)i)W{x + (j + l)i)Z(x - (j + 3)i)W{x + (j + 3)i). 



Comparing the singularity at x = ±(j — l)i and x = ±(j + 3)i we have X ~ Aj-i, Y ~ A 
Z ~ Bj + i and W ~ Bj + \. Then r/j(x) should be factorized as: 



Bj^x) B j+1 (x 
Using the first equation of @ and noting that 



Z 2 
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77,(x) = A,-_i(:e - ji)Aj-i(x + ji)B j+ x(x - (j + 2)i)B j+1 (x + (j + 2)i). (12) 

Considering the singularity at x = ±(j + l)i we have 

Aj^x) _ A j+1 (x) 



(13) 



2-n^x — it x + ie> 
we have 

= exp(^|^ E x _ 2 ^_ ie ) = ex P(^T cot £(* - ^)) > ^o(x) = 1. (14) 
From (0) we have 

A 2j -(x) _ Mx) _ A 2j+1 (x) _ Mx) _ 

B 2j (x) ~ B (x) ~ A ° [Xh B 2j+1 (x) ~ B 1 (x) ~ Ao[X} - (15j 

One can show that A (x) = A' (x). Consider T] 2 - +1 (x + (2j + This quantity approaches 
in the limit of infinite j. Of course near the singularity at 2i, 0, — (4j + 2)i and — (4j + 4)i 
the deviation from becomes big. Nevertheless the region on the complex plane where 

\V2j+i( x + (2j + 1)01 > e is expected to be narrower as j goes to infinity. So we have 

l = Um l+^+i0r + (2j + l)z) = ^ A' {x) A%{x + (4j + 2)i) 



x 



+ (2j + l)i) i-°o A (x) A (x + (4j + 2)z) 
B 2j+1 (x)B2~^(x + (4j + 2)i) B 2j+1 (x - 2i)B 2 ~^{x + (4j + 4)i) 



B 2j (x)B 2j (x + (4j + 2)i) S 2j+2 (x - 22) J B 2j+2 (x + (4j + 4)i) 



(16) 



The first fraction in the big bracket should go to e h / T and the second should go to e h l T . The 
bracket in flT6|) becomes 1 in the limit of infinite j. Then we have 

A' (x) A~ (ioo) 



A (x) A' [ 



a. 



too 



As A (ioo) / A' (ioo) = A (— ioo)/A' (— ioo) = 1/a, we have \a\ = 1. If we choose the phase 
factor a is 1, the ratio of Aj(x) and Bj(x) is always A (x). From (|T0|) and (0) we have 



Aj-i(x - ji)Aj-i(x + ji)B j+1 (x - (j + 2)i)B j+1 (x + (j + 2)i) + 1 

= Aj(x - 3i)A~(x + jQBjix - {j + 2)i)W j (x + {j + 2)i). (17) 

At j = 1 we have 

B t {x - i)Bl{x + i)B a (ac - 3i)~B~ l (x + 3i) = — -= + B 2 {x - 3i)B 2 (x + 3i). (18) 

A {x - t)A (x + i) 

B\, B2 are unknown functions. But this equation and condition 

lim B 1 {x-i)B\{x + i)B l {x-3i)~B~ 1 {x + 3i) = (2 cosh h/T) 2 , (19) 

x— *ioo 

are sufficient to determine these functions. Put 

u{x) = B 1 (x-2i)B[(x + 2i). (20) 

Equation (p~8| ) is written as 

"(* + <)= ,, WtT7 + 3f > . (21) 

v4 (x — zjy4 (a; + i)u[x — zj m(x — 1) 

The l.h.s. has singularity at i, —3i. The first term of r.h.s. has at i, —i, 3i. The second term of 
r.h.s has at 3i, —3i, —i. Assume that u(x) is expanded as follows 

h 00 r- 00 T 

h ^(x- 2nQ - 2%)i fy- 

Consider the contour integral around x = i. Coefficients Cj is determined by 



u(x) = 2 coshQ) + VV ^ + VV 7 ^ -. (22) 

KJ V h Z ?(x-2nQ-2i)J (x - 2nQ + 2i)> 1 ; 



Cj = & {X - l)3 " ^ = / =J^ ^. (23) 

r A (x -i)A (x + i)u(x -i)2ni J A (y)A (y + 2i)u(y) 2ni 



The first sum of r.h.s. of (|22"|) is 

- f exp[-^Ma 1 (y + .)] yJ -i d y _ ^expt-^a^y + i)] 1 dy 



j ./ ., (x — 2nQ — 2%y u(y) 2ni J „ x — y — 2nQ — 2z w(y) 27ri 

. / ^ 4 ( ,_ s _ 2!)exp| _^ aife + !)] _i.^. (24) 

The second sum is calculated in similar way. Thus we get (|5|). From equation (|l0|) we have 

/27rJsinh0 . ,\ . 
1 + T]i{x) = expy — a.i(x))u{x + i)u(x - i). (25) 



Substituting this into eq.(f|) we get eq.(|6|) 

rQ 

f = -T s(x)[\nu(x + i) + \nu(x - i)]dx = -Tlnw(O). (26) 

J — CJ 



3. Analytical solutions 

Ising limit 

In this limit we put A — > oo, J = J z /A. As <fi goes to oo, 2Q becomes 0. Then function 
27r y^ ai(j/) is at \Qy\ > 1 and J z /T at \Qy\ < 1. Function u{y) is 2cosh/i/T at |%| > 2 
and also a constant w(0) at < 2. From eq.® we have 

u(0) = 2 cosh( h/T) + 1 ~ eXP i" J2/T) . (27) 

u(0) 



So we have tt(0) = cosh(/i/T) + ysmh. 2 (h/T) + exp(—J z /T) and known free energy of the Ising 
model. 



Kuniba, Sakai and Suzuki || showed that Takahashi-Suzuki equation for |A| < 1, h = 
can be derived from the quantum transfer matrix and its fusion hierarchy matrices. We are 
confident to derive eqs. (^ ||) from the quantum transfer matrix and its fusion hierarchies. De- 
tails will be published elsewhere. 
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